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j^ Optimum ground states are constructed in two dimensions by using so called ver- 

tex state models. These models are graphical generalizations of the well-known matrix 
product ground states for spin chains. On the hexagonal lattice we obtain a one- 
parametric set of ground states for a five-dimensional manifold of S = § Hamilto- 
nians. Correlation functions within these ground states are calculated using Monte- 
Carlo simulations. In contrast to the one-dimensional situation, these states exhibit 
a parameter-induced second order phase transition. In the disordered phase, two-spin 
correlations decay exponentially, but in the Neel ordered phase alternating long-range 
^. , correlations are dominant. We also show that ground state properties can be obtained 

00 ■ from the exact solution of a corresponding free-fermion model for most values of the 

r**"* , parameter. 
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1 Introduction 

In recent work [1-4], we have presented optimum ground states for quantum 
mechanial spin-1 and spin-| chains with nearest neighbour interaction. Opti- 
mum ground states are exact ground states which are simultaneously ground 
states of the local interaction. They do not involve any excited local states. 
These optimum ground states were constructed using multiple products of fi- 
nite matrices. In |J] we also presented an alternative graphical representation of 
the constructed state, which we called vertex state model. This representation 
has the advantage that it can be easily generalized to higher dimensional spin 
systems. 

In the present work we apply the concept of vertex state models to construct 
optimum ground states for spin-| on the hexagonal lattice. The set of global 
ground states is parametrized by an anisotropy parameter, which determines 
most properties of the ground state. In contrast to the optimum ground states 
for spin chains, this two-dimensional model exhibits a parameter controlled 
second order phase transition. 

The paper is organized as follows: In section @ we present a general parametriza- 
tion of all nearest-neighbour Hamiltonians for £ = I that obey a certain set 
of minimum symmetries. Section S contains a detailed explanation, how non- 
trivial optimum ground states can be constructed for these Hamiltonians using 
vertex state models. The properties of the constructed set of optimum ground 
states are investigated in section H using Monte-Carlo simulations. The most 
important result is the existence of a parameter-induced second order transition 
from a disorded to a Neel ordered phase. In section^ we show that ground state 
properties are asymptotically given by the solution of a free-fermion model. The 
exact solution of this model confirms our numerical results with high precision. 
We finally give a short summary of our results in section pi 

2 Spin- 1 Hamiltonians 

We shall construct optimum ground states for spin-| on the hexagonal lattice 
with periodic boundary conditions. In our previous paper [H we have already 
presented the most general Hamiltonian with nearest neighbour interaction 

ff = £*••„■ (i) 

for S — | systems with the following set of minimum symmetries: 

• Homogeneity: All local Hamiltonians hi_j are equal, they only act on 
different spin pairs. 

• Parity invariance: hi t j commutes with the parity operator Pij, which 
interchanges the two neighbour spins i and j. 



• Rotational invariance in the xy-plane of spin space: hij commutes with 
the pair magnetization operator S* + Si- . 

• Spin flip invariance: S z — > —S z leaves hij unchanged. 

It turns out that all allowed local Hamiltonians can be written as 

hij = A3(|«8)<«8| + |t;_a)<v_3|)+ 

>t{\4){vi\ + \v± 2 ){v± 2 \)+ 
A^(|^)(uj| + bZ 2 )(t;Z 2 |)+ 
*tdKi)(vti\ + \v±n)(v±n\)+ (2) 

Ar(K"X«rl + l«IiX«Iil)+ 

Vl 1 W 01 ) (%L I + A 02 I V 02 ) < W 02 I • 

The local states v? n are defined in Table [l], where we have used the following 
notation for the canonical spin-| basis 

S?\3) = ||3) S?|l) = ||1) 

5f|3) = -§|3) S||T) = -||I). 

In Table [j], /i is the pair magnetization, i.e. the eigenvalue of Sf + S? and p 
is the eigenvalue of the parity operator Pjj . States with p = 1 are symmetric 
under site transposition, the others are antisymmetric. 

For convenience, we fix the local ground state energy at zero leading to a 12- 
parameter model. Because this includes a trivial scale, the number of non-trivial 
interaction parameters is 11. 

A detailed explanation of the concept of an optimum ground state can be 
found in B. This special type of global ground state |$o) satisfies 

hi,j |^o) = for all nearest neighbours i,j , (4) 

i.e. it is also a ground state of all local Hamiltonians. 

3 Vertex state model construction 

Trivial optimum ground states are easily obtained by using simple tensor prod- 
ucts of single spin states \s1S2) for nearest neighbours. Because of parity in- 
variance, this also requires \s2S1) to be a local ground state, which leads to the 
ferromagnetic state of the form 

|*f) = |sisi...si) (5) 

and its antifcrromagnetic counterpart 

\*af) = n \»& n 1*^ * ( 6 ) 
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Table 1: Eigenstatcs of the local Hamiltonian h 
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if the sets A and B define a bipartite decomposition of the lattice and s\ ^ S2- 
In the following, we focus on more complicated ground states. Local ground 
states shall never be simple tensor products of the form \s1S2), instead they 
shall consist of at least two terms 



l^ : 



\S1S2) 



C|S1-1,S 2 + 1) 



(7) 



This special form ensures local S z conservation. In order to be an optimum 
ground state, the global ground state | \&o) has to be annihilated by all local 
Hamiltonians hij. Therefore |^o) contains only local ground states, no local 
excited states are involved. 

The global ground state |^o) is constructed using vertices 




= \si((J,i, 1*2,1*3)) 



(8) 



which map the discrete bond variables i*\ , 1*2, A*3 to a single spin state \s\ (1*1 ,1*2,1*3))- 
The concatenation of two such vertices on the hexagonal lattice is defined by 



taking the tensorial product of the single spin states and summing over the 
shared bond 




(9) 



= E M l S l(Ml,M2,M)>« 

= |si 2 (m.M2,M3,M4)) 



|S2(/X,A*3,M4)) 



The result is a two-spin state which depends on the bond variables /ii, . . . ,^4. 
Because this concatenation can be done for an arbitrary cluster of sites we can 
use it to construct a global state l^o) on the complete hexagonal lattice. I^o) 
does not depend on any bond variables since periodic boundary conditions are 
imposed and thus all /x, are summed over. This construction is much like a 
standard vertex model, except that the 'weights' of the vertices are single-spin 
states instead of numbers and the generic product of numbers is replaced by the 
tensor product of spin states. We therefore call it a vertex state model. 

In order to be an optimum ground state, the constructed global state \^o) 
has to be annihilated by all local Hamiltonians hij (m. Since hij only acts on 
the lattice sites i and j, it is certainly sufficient to demand 



%j 




= 



(10) 



for all fix, . . . ,/i4 and all neighbour sites i,j. In words, the local Hamiltonian 
should annihilate all two-spin states, that are generated by vertex pairs. This 
will ensure that only local ground states enter the global state l^o)- 

Since we are interested in the ground state of an homogeneous global Hamil- 
tonian, we restrict ourselves to vertices which do not depend on the position in 
the lattice. A vertex will be uniquely defined by the state of its bond variables 
and the sublattice it belongs to. 



A particularly interesting optimum ground state is generated by the following 
set of vertices 








X = ID 
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(11) 



In this case, the bond variables \Xi can only take 2 different values, which are 
denoted by arrows on the bonds, a is a continuous real parameter, but a can 



only take the values ±1. These vertices generate the following 9 two-spin states 
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which have to be annihilated by hij. It turns out that this can be achieved by 
putting the following constraints on the parameters defined in (El) 



^2 — A n — ^i — ^oi — ^02 — \i — 
b = c = aa 2 
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(13) 



^3) ^2 

a arbitrary, a = ±1 



Under these conditions, the Hamiltonian still contains 5 continuous and 1 dis- 
crete parameter (c), but the constructed state only depends on a and a. We 
like to point out that the local degeneracy, i.e. the number of local ground 
states, is 9, but the global ground state is unique in the sense that it is the only 
ground state of the global Hamiltonian (fy) for any finite number of lattice sites. 
Appendix A contains a rigorous proof of uniqueness on the finite hexagonal lat- 
tice. This means that condition ( |l0| ) is not only sufficient to obtain an optimum 
ground state, but is also necessary. 



4 Properties of the state 

4.1 General properties 

By construction, the global ground state is an eigenstate of S* otal . In addition, 
it is invariant under a global spin flip S z — ► —S z . Therefore, the single spin 
magnetization (S z ) as well as the global magnetization (Sj: otal ) vanish, which 
indicates an antiferromagnet. 

In the special limiting case a — ► oo, only 4 of the 16 vertices contribute to 

1 Common factors are divided out. 



the global state: 

era 1 3) )E <C — <ra|3) 





(14) 
oa\:-l #^- : |3) . 



With these 4 vertices there are only 2 possibilities of filling the complete lattice: 
A |3) on the first sublattice and a |3) on the second one and vice versa. The 
global ground state is therefore a superposition of both possible Neel statesg 

| Neel! ) + a N ' 2 | Neel 2 ) . (15) 

In the next subsection we present numerical investigations which indicate a Neel 
type behaviour not only in the limit a — > oo but also for finite values of a. 

Another very important special case is a = —y/3, a = — 1. At this isotropic 
point the constructed optimum ground state coincides with the so called valence 
bond solid (VBS) ground state presented in H, rc|. The local eigenstates are 
simultaneously eigenstates of (S^ + Sj) 2 , so the local Hamiltonian can be chosen 
to have complete SO(3) symmetry 
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hi,j = S, ■ S, + — (S< • S,) 2 + — (S< • S,) 3 , (16) 



which is simply the projector onto all states with (S^ + Sj) — 3(3 + 1). The 
authors show in Q] that the global state has exponentially decaying correlation 
functions and conjecture the existence of an energy gap. These results are 
consistent with our numerical investigations. 

4.2 Transformation to a classical vertex model 

In order to obtain correlation functions and other expectation values of the 
constructed optimum ground state, we have to investigate the inner product, 
i.e. the square of its norm. Sticking to the graphical language, the inner product 
consists of two hexagonal lattices placed on top of each other, one for the bra- 
and one for the ket-vector. As a vertex at a specified site only generates single- 
spin states for this particular lattice site, only vertex pairs at corresponding 
sites]] contribute to the inner product. So we just have to take the inner product 
at corresponding lattice sites. The resulting model is a classical vertex modeln 
with 6 emanating bonds as shown in Figure pi The inner product of the global 
ground state is obtained by summing over all bond configurations of this classical 
vertex model. 

2 The overall factor a N is divided out. 

3 i.e. those that are directly on top of each other 

4 i.e. the vertex weights are numbers 



As the inner product of two different S z eigenstates is zero, only 20 of all 
64 classical vertices have a non- vanishing weight. The rule that selects these 20 
vertices is the following: The number of outgoing arrows for the bra- and the 
ket-vector lattice must be equal. There are two special vertices in this classical 
20- vertex model. The vertex with all 6 arrows pointing outwards and that with 
all 6 arrows pointing inwards carry a weight of a 2 . These vertices are generated 
by the inner product of two S z — | or two S z — — § vertices, respectively. All 
other non- vanishing vertices carry a weight of 1. So in the special case a = 1 
the weight of all non- vanishing vertices is 1, which corresponds to an infinite 
temperature in the classical model with the usual interpretation of a weight as 
'Boltzmann weight'. We shall see later that this point lies deeply inside the 
disordered phase of the model. 

4.3 Numerical results 

Even for special values of the parameter a no exact solution of the classical 
vertex model is known. So we have to fall back on numerical methods in order 
to calculate properties of the constructed global ground state. As all vertex 
weights are non-negative, they can be interpreted as Boltzmann weights of a 
thermodynamical system which can be handled using a Monte-Carlo- Algorithm 
(MCA). According to our simulations the system exhibits a phase transition at 
a 2 . = 6.46 ±0.03. For a < a c the S^S^-correlation function decays exponentially 
with a correlation length of order 1. Hence the global ground state is disordered. 
This phase corresponds to the Haldane phase of one-dimensional integer spin 
chains. Both the isotropic point (a 2 = 3) and the infinite temperature point 
(a = 1) lie in this phase. But in the case of a > a c the system is dominated 
by the Neel contributions (flq) and thus the longitudinal correlation function 
shows an alternating long-range behaviour. Figure @ shows a typical correlation 
function for each phase. We have restricted the plot to only one sublattice in 
order to suppress the effect of the alternating sign. 

To determine the type of phase transition, we have investigated the prob- 
ability distribution of Boltzmann weights directly at a c . In our case, there is 
only one peak in this distribution, which clearly indicates a second order phase 
transition. For a first order transition, the phase coexistence would yield at 
least two peaks in the probability distribution M. 

Consistently, the longitudinal correlation function decays algebraically at the 
critical point. Fig. shows the measured correlation function for a system with 
2 x 60 x 60 sites. Again, we have restricted ourselves to only one sublattice to 
avoid the alternation effect. The figure also shows the function 

which connects the measured points nearly perfectly if cj is fitted to 1.62 ±0.01. 
The second term in ( |l7| ) symmetrizes the function, which compensates the effect 
of periodic boundary conditions in first order. 



5 Asymptotic equivalence to a free-fermion model 

5.1 Diagonal vertices and transformation 

Instead of investigating further critical properties using numerical data, we will 
now show that the model is asymptotically equivalent to an exactly solvable 
model for most values of the parameter a. 



The classical vertex model introduced in subsection 4.2 consists of two differ- 
ent types of vertices: If all 3 arrows of a vertex in the bra- vector lattice are equal 
to those in the ket-vector lattice, we call this vertex diagonal. Conversely, off- 
diagonal vertices contain at least one unequal arrow paim. On each sublattice, 
there are 8 diagonal and 12 off-diagonal vertices with non-vanishing weights. 
The crucial result of the Monte-Carlo simulations is that the probability of an 
unequal arrow pair on a bond decays exponentially with increasing values of 
a 2 . As shown in Fig. 0, pj_j is already smaller than 10~ 4 for a 2 = 5. Since 
a 2 > 6, off-diagonal vertices are completely negligible in the regime of the phase 
transition. Properties exactly at the critical point should be independent of the 
model details anyway. This leads us to the diagonal vertex model, which is of 
course a model with only one arrow on each bond. This model on the hexagonal 
lattice consists of the vertices 





(18) 



and identical vertices on the other sublattice. We shall now show that this 
model is a free-fermion model and thus exactly solvable. 

The first step is to reinterpret the hexagonal lattice as a square lattice with 
an elementary cell that consists of two neighbour vertices. Such a pair of ver- 
tices can be combined to a single one by simply summing out the connecting 
bond. The result is a 16-vertex model on the square lattice with the following 
properties: 

5 Actually, these vertices contain exactly two unequal arrow pairs, otherwise the corre- 
sponding vertex weight is zero. 



1. All 16 vertices have a non-vanishing weight. 

2. It is a zero-field vertex model, i.e. each vertex is invariant under a flip of 
all arrows. 

The second property enables us to transform the model to an 8-vertex model 
using the following orthogonal transformation. 

Let us denote the vertex weights of the 16-vertex model by w(ii,i2,i3,ii) 
in which the i v can take the value 1 for a left/down arrow or 2 for a right/up 
arrow. We now define a new vertex model on the square lattice 

w'(jlj2,J3,J4)= (1Q , 

using the orthogonal 2 x 2-matrix 

(20) 




If we concatenate two w'-vertices, both u-matrices on the connecting bond can- 
cel out, because u is orthogonal. The w'-vertex model is therefore completely 
equivalent to the original ui-vertex model. Due to the zero-field symmetry of 
the ^-vertices, only those w'-vertices with an even number of outgoing arrows 
have a non-vanishing weight, i.e. we are left with an 8-vertex model. If the 
usual convention for vertex numbering in the 8-vertex model is used (see e.g. 
[|S| or |9|), the weights of the transformed model can be summarized as 

Wl = |(3 + a 2 ) 2 

w 2 = w 3 = w 4 = w 5 = w 6 = |(-l + a 2 ) 2 (21) 

c 7 =c 8 = ±(l + a 2 ) 2 -2. 

Fortunately, these weights fulfil the free-fermion condition 

UO1U2 + W3W4 = uj 5 uj 6 + uj 7 lu s (22) 

for all values of a, which means the model is exactly solvable. 

5.2 Solution of the free-fermion model 

A detailed investigation and solution of free-fermion models has been done by 
Fan and Wu in [|9| using the method of dimers. Applying their general solution 
to our specific model, the partition function is given by 

1 /•2t 

\nZ = —J h>.{A(<p) + y/Q(y)) dtp, (23) 
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in which 

A(<p) = 

±(a 8 + 18a 4 + 24a 2 + 21) + 2(a 4 - l)(a 2 - 1)cosvj 

Q(<p) = (24) 

(2(a 4 - l)(a 2 - 1) cos<^ - ±((a 2 + l) 2 + 4) 2 ) 2 - 
2( a 4 + 3)(a 2 + l)((a 2 + l) 2 -4) 2 . 

We wish to consider two special points. The first one is of course a 2 = 1. As 
stated earlier, this corresponds to infinite temperature, because all weights of 
the original iu- model are equal to 1, i.e. disorder is maximized. In this context, 
we have to clearify the role of the orthogonal transformation u as a kind of 
duality transformation. At a 2 = 1 only uo\ out of the vertex weights (|2lj) of 
the u/-model is non- vanishing, so in the language of the to'-model, temperature 
is exactly zero. Conversely, the limit a 2 — > oo corresponds to T = in the 
original w-model and to T — > oo in the w'-model. So the u transformation 
exchanges high- and low-temperature limits. Of course, the structure of the 
physical quantum spin ground state is only reflected by the original model, since 
the single spin state at a given site is determined by the number of outgoing 
arrows of the corresponding w- vertex. 

The second important point is defined by Q(<p) = 0, where the system ex- 
hibits a second order phase transition as explained in M . Solving this condition 
for a 2 yields 



= 3 + V12 = 6.464... (25) 



which is in perfect agreement with our numerical result. 

A note added in proof is contained in || , which states that this free-fermion 
model is actually equivalent to the Ising model on an anisotropic triangular 
lattice. This equivalence is capable of explaining the exponent r\ = | in the 
decay of the critical correlation function (|17]). Generally, the correspondence 
between arrow variables a, of the vertex model and Ising spin variables ay is 

a* = ay oy + i , (26) 

in which i' and i' + 1 are the sites on the Ising lattice that 'enclose' the vertex 
model bond i. We can now express correlation functions of the vertex model in 
terms of Ising model correlations 

(oo Or) = (ao o-i u r i ay + i) . (27) 

The leading contribution to this 4-spin correlation function should be 

(o- a-i oy ov+i) ~ (cr Q a r } 2 ex r~= , (28) 

since the critical 2-spin correlations of the Ising model decay asr - *. This result 
is completely consistent with equation (p7|), which is a surprisingly good fit to 
the simulated critical correlation function. 
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6 Conclusion 

We have constructed a one-parametric set of optimum ground states for a five- 
dimensional manifold of spin- 1 Hamiltonians with nearest-neighbour interaction 
on the hexagonal lattice. These ground states arc antifcrromagnetic in the sense 
that the global magnetization vanishes for all values of the continuous parameter 
a. Uniqueness of the ground states on finite lattices with periodic boundary 
conditions can be proven rigorously using an induction technique. For special 
values of the parameters, this ground state coincides with the VBS ground state 
introduced in fl|, m. 

In order to investigate the properties of the constructed ground state we 
have implemented a Monte-Carlo algorithm. The result of the simulations is 
that there is a parameter-induced phase transition from a disordered phase 
with exponentially decaying correlation functions to a Neel ordered phase with 
long-range correlations. This phase transition at a 2 , ss 6.46 is of second order. 
The longitudinal two-spin correlation function decays algebraically at this point 
with an exponent of r\ — h. 

The norm of the constructed global ground state can be regarded as the 
partition function of a classical vertex model on the hexagonal lattice with 
4 bond states. One result of our numerical simulations is that this 4 bond 
state model can be reduced to a free-fermion model with only 2 bond states in 
good approximation for a 2 > 2. The exact solution of the free-fermion model 
reproduces our numerical results with high accuracy. 

Appendix A: Proof of uniqueness 

To prove the uniqueness of the optimum ground state that we constructed in 
section on the finite lattice, we start with an 'empty' hexagonal lattice with no 
spins at the sites. Now we gradually fill this lattice by adding one vertex after 
another (as defined in (O)) at arbitrary empty sites. Whenever two neighbour- 
ing sites are occupied (periodic boundary conditions are imposed), we sum over 
the arrow configurations of the bond between them. After having added several 
vertices and the lattice is partially filled, we have two types of bonds: Internal 
bonds between two occupied sites, where we have summed over, and external 
bonds between an occupied and an unoccupied site with a variable arrow on 
it. During the process of filling, the lattice might exhibit arbitrarily shaped 
and possibly unconnected clusters. Finally, after occupying all sites, we have 
summed over all bonds and end up with exactly the desired global optimum 
ground state. 

This step by step construction enables us to prove the uniqueness of the 
global ground state by induction according to system size. After each step, 
i.e. every time we have added one vertex, we prove the following induction 
hypothesis: 

Running through all arrow configurations on the external bonds 
yields a set of global ground states, which generates all global ground 
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states of the current spin distribution. 

Linear dependencies are only introduced by topologically equivalent 

external bonds.pl 

Since there are no external bonds on the completely filled lattice, we shall have 
proven uniqueness of the constructed optimum ground state. 

For the very first step it is immediately clear from the list of vertices ([II]) 
that the induction hypothesis holds, because we can realize all single-spin states 
on both sublattices. 

Now we assume that the hypothesis holds for systems of size N (and all 
smaller sizes) and add a single vertex. We have to consider 4 different cases: 

1. The new vertex is isolated, i.e. it has no occupied neighbours. 

In this case the hypothesis is trivial. There is no local Hamiltonian that 
connects the new spin with any other spin. So the ground states of the 
iV+1 spin system are simply all tensorial products of N spin ground states 
and all single-spin states. Because we can obviously generate all single- 
spin states with an isolated vertex, we are immediately done by using the 
assumption for N spins. Degeneracy is also trivial. 

2. The new vertex has exactly one occupied neighbour. 

Let «o denote the site where the new vertex has been added and i\ the 
site of its neighbour. Now a new local Hamiltonian hi ^ 1 comes into play, 
which has to annihilate the global ground state. In order to prove the 
induction hypothesis for the N + 1 spin system, we split it into 3 parts: 
The vertex \ip)f at site i , the vertex |<^)f 1:7 at site z'i, and the remaining 
vertex state model for N — 1 spins |$) /32 -' y '. In this notation, the ground 
state condition is 

hi 0til £ ^M^f-W-^O (A.l) 

a,/3i,/32,7 

The meaning of the indices is the following, a represents all bonds of the 
new vertex at iq. Pi are the external bonds of the vertex at i\ before the 
new vertex has been added. 02 are the external bonds of the remaining 
vertex state model. 7 represents the internal bonds between i\ and the N— 
1 spin system. A^ 1 ^ 2 is a tensor of numbers which has to be determined. 

Because hi 0t i x only acts on i$ and i\, the idea is to eliminate |$)^ 2 > 7 from 



(A.l). Normally, we can only do that if the states |<I>)^ 2 ' 7 constitute a basis 
of N — 1 spin space, in particular, they have to be linearly independent. 
In our case, we can still drop |<I>)' 32 ' 7 , because of the following reason: 

The linear dependencies can be eliminated just by restricting the values 
of the $2 index, but leaving the 7 bonds untouched. Figure g shows all 



6 We call external bonds topologically equivalent if they are emanated by the same lattice 
site. A pair of unequal arrows on 2 topologically equivalent external bonds can be swapped 
without changing the global state. 
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neighbours of a new added vertex. These neighbours are all linked to 6 
different sites. So each vertex never emanates more than one 7-b0nd.Fl 
Since linear dependencies can only enter through topologically equivalent 
external bonds (induction hypothesis), only the fa bonds have to be re- 
stricted. We can therefore exclude all over-counted states by restricting 
the A^ 1 '^-tensor, which means no loss in generality, because we are only 
interested in states, not representations. 

After dropping the |$)^ 2 ' 7 , we are left with the following condition 

Km E ^KK 1 ' 7 = °- fOT a11 &>7 , (A.2) 

which is nothing but the ground state condition for the two-spin system. 
Of course, this condition yields exactly the two-spin states ( [L2| ) as its 
solutions, because the local Hamiltonian hi 0i i t is constructed to annihilate 
exactly these 9 local states. As our vertex state model generates these 9 
two-spin states, the first part of the induction hypothesis is proven. 

To check the statement about linear dependency we decompose the system 
into two parts: The two-spin system i$, i\ and the remaining N — 1 spin 
system. We have already proven that the interface between these two 
parts - the 7 bonds - cannot introduce linear dependencies. So the only 
remaining sources are the external bonds of the whole N + 1 spin system, 
which immediately yields the induction hypothesis. 

The above considerations are only valid for o^ and a ^ 00, otherwise 



equation ( |A.2[ ) has more solutions than two-spin states are generated by 
vertex pairs. Consequently the global ground state is highly degenerate if 
a = or a — » 00. 

3. The new vertex has two occupied neighbours. 

In this case, the situation is only slightly more complicated than in the 
previous case. Instead of only one neighbour i\ we now have i\ and i%. 
Consequently, h^^ and hi Q ^ 2 are the new local Hamiltonians that have 
to annihilate the global stater]. The ground state condition is 

(h. . -i-h. . 'i V A0i,02,0 3 . 

Vh Q ,n T ih ,i 2 ) Z-,a,/3i,/3 2 , ,83.71,72 <* (A 31 

l^ ) b>f 1 im |^>f 2 2 ' 72 | $ > /33 ' 7 " 72 = - 

Since the argument of the previous case holds for any number of neigh- 
bours, we can eliminate |$)P3>7i>72 anc j are i e ft w jth 



(h. . + h -IV A0x,02,03 . 

l<l< 1,7 >>f 2 2 ' 72 = 0, for all fa, 71,72 , 



Pi,-n\,„\02,72 _ « <■ " ^ \ ■ ) 



7 This is always true on a bipartite lattice. 

8 «l and «2 can never be crosslinked, because they belong to the same sublattice. 
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which is the ground state condition for a 3-spin system. But this system 
can be constructed using only case 1, so we have already proven that we 
can generate all ground states of this 3-spin system using our vertex state 
model. 

Balance of linear dependencies is checked exactly the same way as in case 
1. Our induction hypothesis is therefore also valid in the case of two 
occupied neighbours. 

4. The new vertex has only occupied neighbours. 

Because the argument is completely analogous to case 2, we omit the 
corresponding ground state conditions here. We end up with a 4-spin 
system, which can be handled completely by using case 1. 

We are now able to generate all ground states for every spin distribution on the 
finite hexagonal lattice with periodic boundary conditions by running through 
all arrow configurations on the external bonds. Because the completely occupied 
lattice does not contain any external bonds, we have proven uniqueness of the 
global ground state that we constructed in section || on the finite hexagonal 
lattice. 
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Figure 1 : Two vertices at corresponding sites combine to a classical vertex with 
6 bonds via the inner product in spin space. 



Figure 2: Longitudinal correlation function in the disordered (dashed curve) 
and in the Neel phase (full curve). 



Figure 3: Longitudinal correlation function at the critical point (dots) and the 
fitted algebraic function (UT 



Figure 4: Probability of unequal arrow pairs as a function of a 2 for a system of 
24 x 24 sites (full curve) and 60 x 60 sites (dashed curve). 



Figure 5: A newly added vertex (unfilled circle) and its neighbours. Note that 
all neighbours are linked to different sites of the remaining lattice. So the bonds 
between the neighbours and the remaining lattice (thick lines) can never cause 
linear dependencies. 
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